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We investigate a one-dimenisonal Hamiltonian system tliat describes a system of particles in- 
teracting through short-range repulsive potentials. Depending on the particle mean energy, e, the 
system demonstrates a spectrum of kinetic regimes, characterized by their transport properties rang- 
ing from ballistic motion to localized oscillations through anomalous diffusion regimes. We etsablish 
relationships between the observed kinetic regimes and the "thermodynamic" states of the system. 
The nature of heat conduction in the proposed model is discussed. 
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The question of how statistical laws emerge from mi- 
croscopic dynamics has been a subject of interest for a 
long time [1] . Studies of relatively simple dynamical sys- 
tems provide a link between deterministic dynamics and 
diffusion phenomena [2,3]. In particular, a number of 
recent dynamical models aim at understanding the nec- 
essary and sufficient conditions for a system to obey the 
Fourier heat conduction law [4]. These examples cover 
only a part of general problem of how kinetic and thermo- 
dynamic properties emerge from dynamics at the atomic 
scale. 

In this Letter we propose a many-particle Hamilto- 
nian model which exhibits a wide range of mass and 
energy transport regimes, and allows to establish rela- 
tionships between microscopic properties and thermody- 
namical and kinetic phenomena. We demonstrate that 
kinetic properties of the system are determined by one 
controlling parameter, which is a mean energy per a par- 
ticle, e. Varying e one can cover the whole spectrum of 
diffusion regimes, from ballistic motion to "frozen" states 
through anomalous diffusion regimes. We show that the 
observed kinetic regimes are strongly related to the "ther- 
modynamic" states of the system which change from a 
"solid" to "gas" phase as the parameter e increases. 

The model describes a system of N classical particles 
each characterized by coordinate Xi and the conjugate 
momentum pi. The particles interact through a repulsive 
short range potential according to the following Hamil- 
tonian 
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where ^ is a strength of inter-particle interaction and a 
is a width of repulsive core. The particles are located on 
an interval of length L, and periodic boundary conditions 
are applied. The mean energy and density of particles are 
defined as e = E/N and g = N/L, respectively. 

In the high energy limit, e/A ^ 1, the system in Eq.(l) 
behaves as a gas of freely flying particles, slightly per- 
turbed by weak interactions. In the opposite low energy 
limit, e/A ^ 1, the system forms a crystal lattice with a 
constant ao = In this solid- like state each particle 



oscillates at the minimum of potential well formed by its 
neighbors. In the gas phase the system is uniform while 
in the solid state it has a discrete translational symme- 
try with a lattice constant oq. According to Kolmogorov 
- Arnol'd - Moser theorem, in both limits, the system 
exhibits an ordered dynamical behavior in a region of 
positive measure in phase space [6]. Due to a difference 
in symmetries of the ordered states one can expect that 
a transition between these two states with the change of 
e will occurs through a mixed disordered state. 

To study the dynamics of the system in Eq.(l) we use 
the following algorithm: At the beginning we put all par- 
ticles at an equal distance ao from each other and give 
them a kinetic energy according to the Maxwell distri- 
bution. After that we rescale velocities of the particles 
in order to get a total energy equal to • e. Then, 
for each time step, we integrate the corresponding dy- 
namical equations using a central difference symplectic 
scheme [5]. To extract an information on the structure 
and dynamics we introduce also a sorted array of parti- 
cles, {a::|°'"*(i)}, where particles are renumbered accord- 
ing their actual instantaneous position, j = sort{i). 

In order to illustrate a nature of excitations in the sys- 
tem we show in Fig. la a time evolution of energy dis- 
tribution initiated by a local perturbation at time t = 
that had a form of a kicked group of a few central par- 
ticles. Here particles with high and low instantaneous 
velocities are displayed by light and dark colors respec- 
tively. The light regions are organized into lines which 
correspond to excitations propagating along the system. 
Three types of excitations with different group velocities, 
which are given by an inverse slope of the lines, can be 
distinguished: flying particles, low-energy phonons and 
nonlinear solitonic excitations which are intermediate in 
energy between them. The first dynamical mode, flying 
particles, dominates in the gas state, e/A ^ 1. Nonlin- 
ear excitations and phonons play the main role in the 
opposite solid-state limit, e/A <C 1, (Fig.lc). An energy 
exchange between different modes is most effective in a 
" liquid- like" state, that corresponds to e/A « 1. Inter- 
actions between modes can lead, for example, to a trans- 
formation of flying particle into solitonic excitation (see 
a circle in Fig. la) or to "burning" of flying particle from 
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FIG. 1: (a) Time evolution of the energy distribution initi- 
ated by kicking of few central particles, for the system in Eq. 
(1) (iV = 128, e/A = 0.25). The circle indicates a collapse 
of flying particle into a nonlinear excitation, (b) The same 
system after time t = 5000. The circle indicates a birth of 
two flying particles from the excitation sea. (c) The set of 
propagating solitons for the system in Eq. (1) for e/A = 1. 

a sea of phonons and nonlinear excitations (see a circle 
in Fig. lb) 

In order to describe quantitatively a structure of dif- 
ferent "thermodynamic" states arising in the system, we 
introduce a reduced distance for a sorted set of particles, 
G — ixj°^*{t) — ~ ^o), and calculate the proba- 

bility density function (pdf) ^{£,), — ap < £, < L/2 (see 
Fig. 2). In the low energy limit the pdf $(^) is governed 
by a peak near the point ^ = 0, which corresponds to a 
crystal lattice (a right inset in Fig. 2). In the gas limit 
the density of states scales as (see left inset in Fig. 2) 

$(^) ^ e-'^(«-''o). (2) 

It is natural to assume that in the intermediate, liquid 
state the pdf, $(^), can be presented in the form 

m^'^soUdiO+'^gasiO (3) 

where ^soUdiO (a dashed line in Fig.2) and ^gasiO (a 
thin line in Fig.2) describe a solid-state and gas contri- 
butions correspondingly. The solid and gas contributions 
can be isolated from the net pdf by fitting its asymptotics 
to the exponential distribution, Eq.(2), at negative and 
large positive ^. Then Eq.(3) allows to introduce the 
order parameter, which presents a fraction of solid phase 
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FIG. 2: The pdf <I>(^) calculated for the liquid state, N = 64, 
A — 1, e — 0.25. The thin line shows fitting to $gas(C) in 
Eq.(2). Dashed line corresponds to $soiid{(,) = 1 — ^gasiO- 
Left inset shows pdf ^{(,) in gas state {e/A = 0.1) and the 
rigth inset corresponds to solid state {e/A = 0.1). 

Fig. 3a shows a dependence of the order parameter cf) on 
the energy density e. It exhibits a phase transition from 
the gas to solid state as e/A decreases. The plateau at 
e/A « 1 indicates on the presence of the third, mixed 
state that can be associated with a liquid-like phase. 
Below we consider regimes of mass and energy transfer 
which correspond to the different thermodynamic states. 

Flying particles. Particles can fly only when their ve- 
locities are higher than a threshold value, Vfi ~ \/2Umax, 
which is determined by the height of the effective poten- 
tial created by neighboring particles. The latter can be 
estimated as a barrier in a periodic lattice, formed by 
particles in the solid state. 

Each particle dynamics can be characterized by the 
mean squared displacement {x^{t)), which in the long- 
time limit, follows 

{x\t))^e (5) 

where a = \ for normal diffusion. All processes with 
a 7^ 1 are known as anomalous diffusion [7,8], being sub- 
diffusion for < a < 1 and superdiffusion for 1 < a < 2. 

The numerical results presented in Fig. 3a clearly show 
three different dynamical regimes: (1) Absence of diffu- 
sion for e/A <C 1 (solid phase). In this case particles 
perform small oscillations around the minima of the po- 
tential, and as a result a = 0. (2) A ballistic regime with 
a = 2 for e/A ^ 1 (gas phase). (3) Normal diffusion, 
a = 1, in the liquid- like state where e/A ~ 1. It should 
be noted that the latter regime exists for a finite interval 
of e. For the values of e corresponding to gas-liquid and 
liquid-solid transitions we have found regimes of anoma- 
lous diffusion, namely a superdiffusion in the first case 
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FIG. 3: (a) The dependence of the order parameter, Eq. (4), 
(j> (black circles) and the diffusion exponent a, Eq. (5) (grey 
circles), obtained by fitting of msd, {x'^{t)), for a calculation 
time t = 10*; The histograms for a duration of a single flight, 
*(t), for (b) e/A = 0.2 and (c) e/A = 1.2; iV = 512. 

and subdifFusion the second case (Fig. 3) [9]. 

In the case of superdiffusion each particle performs 
both flights and trapped oscillations, randomly switch- 
ing between these two dynamical modes (see Fig. la and 
b). We have obtained the histogram for a single flight 
time ^'(t), accumulating the lifetime and length dur- 
ing the interval where the sign of the particle veloc- 
ity is fixed. It shows an asymptotic power law decay 
(Fig.Sb), ^(t) - t^T, with the exponent 7 « 2.5 for the 
gas phase {e/A = 10) that corresponds to superdiffusion 
and 7 « 3 for the high energy edge of the liquid phase, 
at e/A = 2, that corresponds to a transition to normal 
diffusion regime. The related pdf has the same power 
behaviour. We note that the flying particles do not have 
a constant velocity. 

The subdiffusion, which has been found in the transi- 
tion region between solid and liquid states, corresponds 
to short random jumps of the particles interrupted by 
long trapping events. At low energies the probability of 
strong fluctuations which can result in a particle escape 
from the potential well becomes very low. This leads to 
anomalously long trapping of particles. The particle mo- 
tion becomes subdiffusive when the pdf of trapping times 
has a power law asymptotic, ^'(i) ^ t^^^f^ with /3 < 1 [8]. 
In this case the diffusion exponent a in Eq. (5) equals 
/3, a = (3 < 1. The subdiffusive dynamics can be quanti- 



tatively described using a continuous-time random walk 
(CTRW) approach [8]. 

Solitary excitations and phonons. In the low-energy 
regime the system tends to build a crystal-type structure. 
Local energy fluctuations can create here only relatively 
smooth nonlinear excitations as well as small-amplitude 
phonons. In the phonon state particles are mainly local- 
ized in the middle point between two nearest neighbors, 
and the variable introduced earlier can be used as a 
small parameter, ^j/ao <C 1 [11]. Lianerization of the 
equations of motion for the Hamiltonian in Eq.(l) gives 

the elastic constant C — 14(1 — '(3^)^ ("6)2^ which 
contains all information on the phonon spectrum [11]. 
This approximation is not valid for nonlinear excitations, 
where particles can strongly deviate from their equilib- 
rium positions. However, in the case of smooth excita- 
tions the displacements of particles with respect to the 
corresponding middle points change slowly with parti- 
cle number j in the sorted array x^"""*, and the differ- 
ence j/j = ix^j'^l -I- x'/°^l)/2 — Xj can be considered 
as a smooth variable y = y{x,t). One can rewrite the 
Hamiltonian in Eq.(l) in terms of the continuous vari- 
able y{x, t) and prove that the solitonic nonlinear excita- 
tions propagating with a constant velocity, y = y{x — vt), 
indeed exist in the system. 

We have found that such solutions exist only for par- 
ticle densities which are below a critical density, Qcrit = 
1/aocrit = v^/ct. For the densities higher than the crit- 
ical one the width of the excitations reduces to a sin- 
gle particle scale and the smooth nonlinear excitations 
cannot exist. Under this condition the liquid state (the 
plateau in Fig. 3a) almost disappears and a transition be- 
tween solid and gas phases goes directly through a " subli- 
mation" process. Numerical studies of the nonlinear solu- 
tions show that they have properties of " quasi-particles" . 
They propagate in both directions practically without a 
loss of energy (see Fig.lc) and collide with each other 
another. Interactions between different types of excita- 
tions produce an intensive energy exchange in the liquid 
state that leads to the most efficient thermalization of 
the system in this state. 

Heat conduction. The proposed system provides a pos- 
sibility to simulate kinetic properties of classical low- 
dimension systems [4]. As an example we consider a 
problem of heat conductivity. For this purpose the sys- 
tem in Eq. (1) is coupled to the heat reservoirs placed 
at the walls x — and x ^ L. The temperatures at the 
left and the right walls are given by Ti and Tr (Ti > Tr), 
respectively. When a particle collides with a wall at tem- 
perature Ti r, it is reflected back with a velocity chosen 
from the distribution f{v) — {\v\/Ti^r)exp{—mv'^/Ti^r)- 
To calculate a temperature profile we evaluate time av- 
erages as follows: we divide the interval L into a set of 
equal cells, X^, s = 1, TV - 1, with a length AAT = L/N. 
In Fig. 4 we show the local temperature distribution 
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FIG. 4: Local temperature profiles Ti for e/A — 10 (open 
circles), e/A = 1 (filled circles), and e/A — 0.1 (trian- 
gles); N = 100, Ti = 2, Tr = 0.1. Inset shows local tem- 
perature profile and local density of particle, calculated for 
Ti = 0.1, = 2, iV = 128, e/A = 0.25. Lines obtained by 
averaging over 100 realizations after a transient time t — 10^. 



T{x) = (Ws)/2 found for different thermodynamic phases 
of the system in Eq. (1). Here (■(;s)/2 is an average kinetic 
energy of the s-cell and ( ) denotes a time averaging, s 
denotes the location along the interval and corresponds 
in a continuum description to x. In the gas Hmit the 
temperature profile exhibits a wide plateau which is typ- 
ical for a ballistic mechanism of conductivity (Fig. 4, open 
circles). In this case the heat transfer is determined by 
flying particles. In the opposite, solid state, limit we also 
find ballistic conductivity which is determined by propa- 
gating quasi-particles: solitonic excitations and phonons 
(Fig. 4, squares) [4]. In the hquid state the T{x) has a 
strong nonlinear profile with a substantial nonzero slope 
that is a result of intensive scattering processes and en- 
ergy exchange between all dynamical modes (Fig. 4, filled 
circles). The local density of particles, p{x) — {ns)/N, 
where Ug is the average number of particles per cell, also 
demonstrates a nonlinear profile, which is complimen- 
tary to the temperature profile. The observed variation 
of p{x) can be explained by an accumulation of slow par- 
ticles at the cold end and by fast escape of flying parti- 
cles from the hot end. A relation between temperature 
and concentration gradients, which is usually assumed 
for a macroscopic system close to equilibrium, naturally 
emerges from the proposed Hamiltonian model. 

Our calculations show that the proposed system can 
be used as a model for thermo- engine that transforms 
heat into a directed motion. For this purpose the en- 
semble of particles should be coupled to a third body, 
i. e. a cargo. Slow particles moving from the cold end 
affect the cargo stronger than fast particles that move 
from the hot end. This asymmetry leads to a directed 
motion of the cargo. Another example of transformation 
of thermal energy into a directed motion is a produc- 
tion of electric current due to temperature differences of 



thermostats. To simulate this phenomenon we have con- 
sidered two kinds of particles with different masses and 
opposite charges. When the difference in the masses of 
the particles is large enough the light ones can fly ballis- 
tically and transfer current while the heavy ones remain 
immobile. A detailed consideration of phenomena related 
to a directed motion will be published elsewhere. 

In summary we have introduced a simple dynamical 
model that establishes the relationships between micro- 
scopic properties of Hamiltonian systems and "macro- 
scopic" thermodynamical and kinetic phenomena. De- 
pending on the energy the model exhibits three well- 
deflned states and a wide spectrum of kinetic phenomena. 
An interesting and still unresolved question concerns a 
scaling of Fourier coefficient for the thermoconductivity, 
k ^ , with a length of the system, L, [4]. The solution 
of this problem requires much longer simulations than 
that reported in this Letter. On the base of our calcu- 
lations we expect that normal heat conduction, /c ~ L'^, 
should be observed in the liquid state. 
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